clear all
ai = 13e-3;
ao = 30e-3;
l = 10e-3;
LorentzForceCase1_Field3D(ai,ao,l);